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Institute of Physics, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia 

We study the implications of the simplicity constraint in the spincnbe model of quantum gravity. 

By relating the edge lengths to the integer areas of triangles, the simplicity constraint imposes 
very strong restrictions between them, ultimately leading to a requirement that all 4-simplices in 
the triangulation must be almost mutually identical. As a surprising and unexpected consequence 
of this property, one can obtain the CDT state sum as a special case of the spincnbe state sum. 

This relationship brings new insight into the long-standing problem of the relationship between the 
spinfoam approach and the CDT approach to quantum gravity. In particular, it turns out that the 
spincnbe model contains properties of both approaches, providing a single unifying framework for 
their analysis and comparison. In addition, the spincube state sum also contains some other special 
cases, very similar but not equivalent to the CDT state sum. 

PACS numbers: 04.60.Pp 


I. INTRODUCTION 

Loop Quantum Gravity (LQG) 0,0 is one of the mainstream approaches to the fundamental problem of the gravi¬ 
tational field quantization. Various formulations of concrete models of LQG can be roughly divided into canonical and 
covariant. The canonical formulation focuses on the Hamiltonian formalism of gravity and its canonical quantization. 
The covariant formulation focuses on providing a rigorous definition to the gravitational path integral 

Z = J 

This is mostly being done by discretizing the integral, like in ordinary quantum mechanics. The models thus obtained 
are generically called state sum models (see Appendix]^, and they describe gravitational and matter fields by suitably 
labeling a piecewise-linear 4-manifold. In other words, the path integral is defined as a sum over certain labels on a 
4-dimensional lattice, most commonly an irregular triangulation. The most developed models of this type are called 
spinfoam models, see for example [2-0 ■ alternative definition of the path integral is to fix the labels to have 
constant values, and perform the sum over some class of manifold triangulations instead. The most developed model 
of this type is called Gausal Dynamical Triangulations approach (GDT), see 0 for a review. As it turns out, the two 
choices of what to sum over in the path integral lead to related but very different models, whose properties are not 
easy to compare. This represents a long-standing open problem of the relationship between the spinfoam approaches 
and the GDT approach. 

Recently, a novel generalization of spinfoam models has aweared, called the spincube model 0. It represents a 
generalization of spinfoams through higher category theory |8|. In particular, whereas spinfoam models are based 
on the structure of a Lie group, the spincube model is based on the structure of a Lie 2-group, a 2-categorical 
generalization of the notion of a group (see 0 for a comprehensive list of references on 2-groups). As a consequence, 
while in spinfoam models one uses the structure of a labeled simplicial 2-complex, in the spincube model one uses the 
structure of a labeled simplicial 3-complex (hence the names “spinfoam” and “spincube”). The spincube model was 
introduced in order to improve certain properties of spinfoam models, like the presence of the tetrad fields, coupling 
of matter fields to gravity and the semiclassical limit, which are issues hard to deal with in spinfoam models. 

It thus comes as a fresh surprise that the spincube model naturally incorporates both the GDT and spinfoam 
approaches within a single model. In this paper we will demonstrate explicitly how and why this happens, thus giving 
some new insight to the relationship between spinfoams and GDT discussed above. This is an unexpected result, 
which came about accidentally, during the study of certain properties of the spincube model. 

In short, what happens is the following. The spincube model consists of the state sum quantization of a topo¬ 
logical BFCG theory 0110, to which one then applies a simplicity constraint. This is in complete analogy with 
spinfoam models, which are constructed by starting from a state sum of the topological BF theory and then applying 
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the simplicity constraint. As it turns out, in the spincube model the simplicity constraint has a simple geometric 
interpretation — the areas of all triangles (as calculated using the edge lengths for each given triangle) must have 
integer values, in certain units. But since there are more triangles than edges in any given triangulation, the simplicity 
constraint is represented by a system of equations which is overcomplete in edge lengths, while taking into account 
integer areas as variables leads to a complicated system of Diophantine equations, which may happen to have no 
solutions at all (this was first noticed in A more detailed analysis of the simplicity constraint equations was thus 
necessary, to study whether they can be imposed at all or not. That analysis is presented in this paper, with a result 
that is doubly surprising. First, despite the complicated Diophantine equations, the set of solutions for the simplicity 
constraint is indeed nonempty, and several classes of solutions will be explicitly constructed. And second, all these 
solutions have a nontrivial resemblance to the CDT model, while one class turns out to transform the spincube state 
sum model precisely into the CDT state sum model. 

One can thus say that the spincube model represents a generalization of both spinfoam and CDT approaches, and 
encapsulates the properties of both into a single model. This is very interesting, as it opens novel avenues for the 
study of both approaches, and especially their mutual comparison. 

The layout of the paper is as follows. In section El we will give a short overview of the spincube model of quantum 
gravity, introduce notation and conventions. An important step in the construction of the model is the way one 
imposes the simplicity constraint. Section IlIIl is devoted to a detailed discussion of the difference between imposing 
the simplicity constraint weakly and strongly. It will turn out that the strong simplicity constraint is nontrivial to 
impose, and might not have any solutions. In section ITVl we give a proof that the strong simplicity constraint actually 
does have solutions, and we construct one class. Another class of solutions is then discussed in section El and we 
show that for these solutions the spincube model reduces precisely to the CDT model of quantum gravity. Section IVTl 
deals with the implications of possible nontrivial spacetime topologies to the spincube model. Our conclusions and 
discussion of the results are given in section Em while the Appendix contains some further technical information. 

The notation and conventions are as follows. The local Lorentz indices are denoted by the Latin letters o, 6 , c,..., 
take values 0,1, 2, 3, and are raised and lowered using the Minkowski metric rjab with signature (—, +, +, +). Spacetime 
indices are denoted by the Greek letters ^,v,..., and are raised and lowered by the spacetime metric g^i, = rjabe^'fj,e^ , 
where are the tetrad fields. 


II. SPINCUBE MODEL 

We begin by giving a review of the spincube model, which was introduced in detail in Q. In this paper we do not 
aim to provide a detailed motivation for the introduction of the spincube model, nor a detailed analysis. Rather, we 
merely repeat the minimum self-contained amount of material given in Q needed to fix the notation and to spell out 
the main properties of the model. A reader interested in more details about the model itself, the Poincare 2-group 
and representation theory of 2 -categories can look up the literature [MS Ed- 

Let AI 4 be a four-dimensional spacetime manifold. The spincube model of quantum gravity is based on the classical 
action 

j + (Bafc-eabcde^Ae''), (1) 

JMi 

which is called the constrained BFCG action for the Poincare 2-group. The first two terms represent the topological 
BFCG action [^, d|, which is a 2-categorical generalization of the standard BF action. The last term is called the 
simplicity constraint. The fundamental variables in this action are two-forms B°‘^, /3“ and one-forms e“, 

The B°‘^, (/)“*’ and e“ are Lagrange multipliers, while uj°'^ and (3°' together constitute the 2-connection for the Poincare 
2-group. They give rise to corresponding two-form and three-form field strengths 

-f a;“e A G“ = d/3“ -f A 

Note that the quantities traditionally denoted Fab and G“ in the BFGG action have been conveniently relabeled 
here as Rab and e“, since they have the interpretation of spacetime curvature two-form and the tetrad one-form, 
respectively. 

The variation of the action © with respect to all of its variables yields the following set of equations of motion: 

: Rab-4>ab = 0, 

: Bab - £abcd A = 0, 

: VRab - e[a A /3b] = 0, 

( 5 / 3 “ : Vea = 0 , 

(5e“ : Ga + 2ea6,d/"Ae^ = 0. 
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If the tetrad field is assumed to be nondegenerate, det(e“;i) 7 ^ 0, this system of equations can be rewritten into an 
equivalent system, consisting of the equations that solve for the Lagrange multipliers, 

4^ah — Rab^ Rab — ^abcd 6 A 6 , 

the equations for /3“ and torsion two-form = Ve“ = (ie“ -I- a;“c A e°, 

pa ^ 

and the equation for the tetrad field e“, 

eabcdR^‘' A = 0, 

which is the differential-form notation for the Einstein equations of general relativity. Thus we see that the action (ED 
is a reformulation of general relativity (or more precisely of the Einstein-Cartan theory, see 0)- 

This reformulation is convenient because it has two key properties. First, the action is split into the topological 
BFCG part and the simplicity constraint. This is in complete analogy with the Plebanski action, which is split 
into the topological BF part and the simplicity constraint. The usefulness of this split lies in the fact that one 
can perform the quantization of the topological sector, and then impose the simplicity constraint, thus defining a 
theory of quantum gravity. In the case of the Plebanski action, this program leads to spinfoam models of quantum 
gravity (see (H)- In the case of the action ED, the same program leads to a 2 -categorical generalization of spinfoam 
models, dubbed spincube model. The difference lies only in the fact that the quantization of BF theory is based on 
the underlying structure of a group (giving rise to a colored 2 -complex, hence “spinfoam”), while the quantization 
of the BFCG theory is based on the underlying structure of a 2-group (giving rise to the colored 3-complex, hence 
“spincube”). 

The second key property is that, in contrast to the Plebanski action, the action ED features tetrad fields as explicit 
variables in the topological sector. After the quantization of the topological sector, the presence of tetrads manifests 
itself as the presence of a new set of colors on the 3-complex, which are interpreted as edge lengths. The spinfoam 
models do not feature these labels, and consequently run into trouble when one tries to couple matter fields, especially 
massive fermions (see Appendix in [l^ for a short overview). In contrast, the coupling of any type of matter fields to 
the spincube model is completely straightforward 0, and represents a clear improvement over the spinfoam models. 

Specifically, the mass terms of matter fields in the discretized action (and also the cosmological constant term) 
depend on the 4-volume of each 4-simplex in the triangulation. This 4-volume term may be hard or even impossible 
to express in terms of the fundamental variables of a given spinfoam model. For example, in area-Regge variables the 
4-simplex is assigned 10 areas for its triangles. These variables do not fix a unique value of the 4-volume (see [l^ and 
also Appendix IHD, so mass terms in the discretized action are not well-defined. Similarly, in the Engle-Pereira-Rovelli- 
Livine/Freidel-Krasnov spinfoam model one uses triangle areas and tetrahedra 3-volumes as variables, which also do 
not fix a unique choice of the 4-volume without the knowledge of at least one edge length 0 , which is not a variable 
of the theory. In the spincube model, this issue is resolved in arguably the simplest possible and geometrically most 
natural way, due to the explicit presence of tetrad fields in the topological part of the action ED- On a triangulation, 
these variables become the edge lengths, so given a set of 10 edge lengths in a 4-simplex, its 4-volume is uniquely 
determined and easily expressed in terms of the Cayley-Menger determinant formula (see (ID2I) 1. Therefore, all mass 
terms in the discretized action are automatically well-defined. 

Let us now summarize the details of the quantum theory corresponding to ED- The spincube model is one particular 
example of a state sum model (see Appendixl^for a short introduction to state sum models), and is defined as follows. 
Let T{M. 4 ) be a triangulation of the spacetime manifold A 44 . Triangulation T contains vertices v, edges e, triangles 
A, tetrahedra r and 4-simphces a. The total numbers of edges and triangles will be denoted as E and F, respectively. 
Each edge e is labeled by a non-negative parameter 4 G R,)", each triangle A is labeled by an integer toa G 2, while 
vertices, tetrahedra and 4-simplices do not carry any additional labels. Thus the set of all rriA represents the 
fundamental variables which describe the gravitational field. Further, the vertex amplitude Ay and the tetrahedron 
amplitude At are chosen to be trivial. 


Av(l,m) = l, Ar(^TO) = I, 

while the 4-simplex amplitude is chosen as 

A — „iST(l,m) 

where So- is the Regge-like action for a given 4-simplex u, 

Sa{l,m) = ^ |TOA|0A,a-(O - 

Ago- 
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Here 0A,cr(O the dihedral angle corresponding to the triangle A of the 4-simplex u, and is a function of the edge 
lengths of tr. As we shall discuss in the next section, the integer rriA is related to the area of the triangle A. 

The state sum of the spincube model is then written as 


Zt= dh... X! ■ " X! n -4e(^TO) TT AA{l,m) TT 

JR+ Jk+ miGZ mpeZeeT 


,iSa{l,m) 


( 2 ) 


AgT 


<tGT 


The product over all 4-simplices can be rewritten in the form 


n< 

<7£T 


,z5a(Z,m) 


= exp 


= exp 


= exp 


= exp 


L <y^T 
i 


i ^ Sa(l,rr. 

creT 

Y l"iA|0A,a(O 


. o-GTAGct 

i Y I"^a| ( Y 

. Aer VcrgA / 


(3) 


i Y^ I™a|(5a(0 

L Aer 


_ piSR{l,m) 


where Ja = So-sa ®a.<t is the deficit angle for the triangle A in Minkowski space, and Sr is the total Regge-like 
action over the triangulation T{Mi). It differs from the usual Regge action by the presence of the integer rriA in 
place of the area of the triangle, Aa{1) — the former is an integer independent of edge lengths, while the latter is a 
function of edge lengths given by the Heron formula. This difference will be discussed in detail in the next section. 
The spincube state sum can thus be written as 

Zt= [^dh... [dlE Y ■■■ Y. (4) 

•^*0 -^*0 miez mpeZeGT AgT 

where the remaining edge and triangle amplitudes Ae and Aa are yet to be specified. Their forms are based on the 
way we choose to impose the simplicity constraint, which is the main topic of the next section. 

One should note that the model can be easily extended to include matter fields, cosmological constant, etc. However, 
in this paper we will omit such extensions in order to make the exposition more transparent. 


III. IMPOSING THE SIMPLICITY CONSTRAINT 

Traditionally, given a classical action, one constructs the state sum model in the following way. First, one rewrites 
the classical action into a topological sector and the simplicity constraint. Second, one constructs the state sum model 
for the topological sector of the theory, so that the resulting state sum is a triangulation independent topological 
invariant. Third, one modifies the amplitudes and imposes the simplicity constraint in a certain way, in order to pass 
from the topological state sum model to the non-topological state sum model which represents the quantization of 
the full classical action. 

The crucial difference between the topological and non-topological state sum models lies in the implementation 
of the simplicity constraint, which represents the difference between the full theory and its topological sector. The 
simplicity constraint is implemented as a suitable relation connecting the independent variables of the topological 
state sum model. In the case of the spincube model, those variables are positive real numbers le and integers ttia, 
associated to edges and triangles in the triangulation, respectively. In the variables of the classical action o the 
simplicity constraint reads 

^ab ^ ^abcd A Cd , 

and is a classical equation of motion obtained by varying the action 0 with respect to (j)ab- Passing to the variables 
le^ruA defined on the triangulation, the simplicity constraint takes the following form: 


'yip I'm A I — Ar (Zjj , j ^£3 ) I 


VA e T{M4) , 


(5) 
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where ei, £ 2 , £3 G A. Here Ah{1) is the Heron formula for the area of the triangle A with edge lengths I, 


Anihih, I 3 ) = \Js{s — li){s — l 2 ){s — I 3 ) , (6) 

where s = (^i + + ^ 3 )/ 2 , while 7 is a dimensionless parameter determining the overall scale, and the unit of length is 

provided by the Planck length Ip. The simplicity constraint has a natural and obvious geometric interpretation — for 
every triangle A in T, the integer ttia determines the area of the triangle, with a universal proportionality constant 

ill- 

We should note that the factor 7 is a free parameter of the theory. It should not be confused with the Barbero- 
Immirzi parameter which is also commonly denoted by the same letter 7 . The Barbero-Immirzi parameter can be 
introduced as a coupling constant if one adds the so-called Holst term [1^ to the action O- In this paper we do not 
introduce it. 

One can immediately see two basic properties of the simplicity constraint. First, it imposes triangle inequalities on 
all edge lengths le in the triangulation. Namely, in triangulations with Lorentzian signature it is perfectly possible for 
the three edge lengths to violate triangle inequalities. This would give rise to an imaginary value of area ®, which 
corresponds to a triangle coplanar with some time-axis in Minkowski spacetime. However, simplicity constraint (O 
requires all triangles to be “spacelike”, i.e. to have real-valued (and positive) areas, which in turn imposes the triangle 
inequalities through dH]). 

Second, the simplicity constraint transforms the Regge-like action from 

SR{l,m) = I"^a|i5a(0 

AgT 


into the proper Regge action 


^Regge(0 = 4 E^a(0<5a(0- (7) 

AeT 

However, there is one big issue with the simplicity constraint ®. Namely, as it was first noted in 0 , the system 
of equations ® is not guaranteed to have any solutions. This can be seen as follows. If the total number of edges in 
T is E, and the total number of triangles is F, then for every triangulation we have F > E, while the equality holds 
only for a single 4-simplex. This means that we have in total F equations for F integer variables and E real variables. 
If we write the system generically as 


|toi| 

— ^ 1 (^ 1 , • 


\mp\ 

= Ap{li,. 


tob+iI 

= Appi{li,... ,Ie 

\mp\ 

= Ap{li,. 

■■Je), 


where Ai,..., Ap are suitable functions, we can in principle solve the first E equations for as functions of mi ,... ^mp- 
Substituting those expressions into the remaining F — E equations, the latter can be written in the form 

|m_E+i| = /i(mi,... ,mE), 


|m_F| = fp^pimi,... ,mE), 


( 8 ) 


where /i,..., fp-p are implicitly defined functions, too complicated to be expressible in terms of elementary functions. 
Given that all mA are integers, equations ® are thus very complicated Diophantine equations, and they are not 
guaranteed to have any solutions at all. Therefore, short of providing the proof that the set of solutions for ([ 8 ]) is 
never empty, the simplicity constraint system may be impossible to enforce as it stands. 

In m an alternative strategy was proposed — to implement the constraint “weakly”, such that it holds only in 
the classical limit of the theory. In particular, it should be written as 


I TO A 


^^2 ^ ) 


VA G TiMi) 


(9) 
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where again ei, £ 2 , £3 G A. Here [... J represents the “floor” function, which returns the integer part of its argument. 
This definition of the simplicity constraint completely circumvents the issue of Diophantine equations. Moreover, as 
argued in [lll |. in the classical limit (defined as an asymptotic expansion when mA,le —^ 00 ) the difference between 
([S]) and ([H]) is always a lower order correction, making the two systems of equations asymptotically equivalent. 

In the formalism of the state sum model, the implementation of the weak simplicity constraint ([9]) is given by the 
following choice of the edge and triangle amplitudes, 


A,{l,m) = l, -4a(?,to) = X ( |toa|- 


Ml 




where xM) i® indicator function, having value 1 if a: = 0 and is zero otherwise. Substituting this into the state 
sum o, one can immediately perform the summations over all toaj and the state sum reduces to the one defining 
the Regge quantum gravity model: 


'T’weak 

^ rp 


dh... dlE , 


( 10 ) 


where the domain of integration over edge lengths is a complicated subset of due to the triangle inequalities 

imposed by the simplicity constraint. The action in the exponent is 


Sr{1) = 


E 

Aer 


1% 


(4ga) 


^a(0 , 


( 11 ) 


which is asymptotically equal to the proper Regge action © in the classical limit —>■ fc —>■ 00 (see Appendix ICl for 

proof). Note that the choice of the edge and triangle amplitudes defines the measure of the discretized path integral, 
and the product of these amplitudes over all edges and triangles can be identified with an appropriate Jacobian 
determinant. 

Nevertheless, one can argue that imposing the simplicity constraint weakly might be unsatisfactory on various 
grounds, and it is legitimate to ask if the constraint can be imposed strongly. Under the assumption that the 
Diophantine system discussed above has at least one solution for l^, denoted /g = Lf^{fh), where the fh denote the set 
of variables unconstrained by the equations, one can implement the strong simplicity constraint ([5]) with the following 
choice of edge and triangle amplitudes: 


Ae{l, m) = S {k - L^{m)) , Aa(^ m) = x ^|wa| - iLeeA{m))^ . 

Substituting into (|T]), we can now integrate over all edge lengths, and evaluate all sums over m except those over rh. 
The resulting state sum is 


Zrp ^ - E giSR(L(m)) ^ ^22) 

fm} 

where the sum is taken over all independent sets of integers rh. The action in the exponent is 

SR[L{m)) = ^ E AH{L^^A{m))M{L{m)) > 

agt 

and it is equal to the proper Regge action © evaluated on the simplicity constraint solution 4 = L^(m). 

The implementation of the strong simplicity constraint rests upon the assumption that the complicated Diophantine 
equations © have a nonempty set of solutions, which is not obvious. It may thus come as a fresh surprise that a 
few classes of solutions can indeed be found. Moreover, the resulting spincube model can then be naturally related to 
a completely independent approach to quantum gravity, namely the Causal Dynamical Triangulations approach Q. 
This is both completely unexpected and a very interesting result. 

In the next section, we will perform an explicit construction of one class of exact solutions of the strongly imposed 
simplicity constraint, and after that we will discuss the relation between the spincube model and the CDT approach 
to quantum gravity. 
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IV. SOLUTION OF THE SIMPLICITY CONSTRAINT 


Let us turn to a constructive proof that the simplicity constraint ([5]) always has at least one solution. The proof 
will be done in three steps. First we will discuss the case of a single 4-simplex, then the case of two 4-simplices sharing 
a common tetrahedron, and finally the general case of an arbitrary triangulation. 

A single 4-simplex has 10 edges and 10 triangles. In that case, the simplicity constraint ([5]) has the general form: 

|wi| = ... ,lio), 

: (13) 

l™iol = —p Aio(Zi,...,/lo) • 

This system of equations always has solutions. A simplest example is an equilateral 4-simplex, with areas and edges 
given as 


mi = • • • = TOio = k , li 


- /iQ — 2/. 



for any fc S Z. For more complicated choices of m’s, numerical analysis suggests that several solutions for /’s may 
exist, since there are ten polynomial equations to be solved. However, note that we require the edge lengths to be 
real-valued and positive, and moreover they must satisfy triangle inequalities. The choice of m’s might be such that 
this is impossible to satisfy, in which case the system does not have any solutions. Therefore, it is important to stress 
that we are not attempting to solve the system for ten Vs given any arbitrary m’s, but rather to solve the system 
for Vs and for m’s simultaneously, within their respective domains. The equilateral example proves that the set of 
solutions is nonempty, and numerical analysis (of Monte-Carlo type) shows that the set of solutions is actually quite 
rich (see Appendix [Pl for details). 

Next we move to a less trivial case of two 4-simplices sharing a common tetrahedron. The 4-dimensional figure is 
depicted on the following diagram: 


1 



The first 4-simplex, cri, is determined by the vertices (I, 2, 3,4, 5), while the second, CT 2 , is determined by the vertices 
(1,2,3,4,6). They share the common tetrahedron r determined by the vertices (1, 2, 3,4), depicted with thick edges. 
There are four triangles of r, shared by both 4-simplices, namely 


(1,2,3), 

In addition, cji contains six more triangles 

(1,2,4), 

(1,3,4), 

(2,3,4). 


(5,1,2), (5,1,3), (5,1,4), 

while (72 contains its own additional six triangles 

(5,2,3), 

(5,2,4), 

(5,3,4), 

(6,1,2), (6,1,3), 

(6,1,4), 

(6,2,3), 

(6,2,4), 

(6,3,4). 


In total, the figure has E = \A edges and F = 16 triangles. For each triangle we write a simplicity constraint equation, 
giving rise to 16 equations: 

|to123| = 


|W634| = 


ill 


A_f/(Zi2, ^13, hs), 


ill 


1^h{163, Iqa, hi) ■ 


( 14 ) 




As discussed in the previous case, one can always solve the system of the first ten equations, for the 4-simplex ai. 
Denote this solution as 


mi23 

= mi23 , 

^12 — Li2 

m534 

= m534 , 

154 = L 54 


(15) 


One could now arguably also solve four of the remaining six equations, for example for triangles (6,1,2), (6,1,3), 
( 6 ,1,4) and ( 6 , 2,3), giving rise to additional four m’s and four L’s. Substituting all previous results into the remaining 
two equations, one gets 


|w 624| =/(m, L), lme 34 l = g(m,L). (16) 

Here / and g are two functions implicitly defined by the substitution of 14 m’s and 14 L’s into the remaining two 
equations. However, as discussed previously, there is a priori no guarantee that there will exist a choice for (m, L) 
such that m 624 and m 634 are integers. But there is a beautiful geometrical argument which proves that this choice 
indeed exists. Namely, given the generic solution (TTSl) for cti, choose the following for (T 2 : 

^61 = A 51 , Iq 2 = A 52 , ^63 = A 53 , Iq 4 , = L 54 . (17) 

In other words, we choose the edges of tT 2 to be equal to corresponding edges of (Ti, making the two 4-simplices identical 
(up to reflection symmetry). Since the triangle areas of tri are integers by construction, so will be the triangle areas 
of ( 72 . In particular, we have 


^612 = ni512 , ni623 = ^523 > 

m6l3 = fh 5 i 3 , m624 = 1^524 , (18) 

m6l4 = fh5i4 , m634 = m534 . 

In this way, all 16 simplicity constraint equations are simultaneously satisfied, giving us one explicit solution for the 
case of two connected 4-simplices. Numerical search for other solutions of the two equations (HU) has proved fruitless, 
suggesting that this is the only possible generic solution. In special cases (like isosceles 4-simplices) there may be 
more solutions, and examples of these will be given in the next section. 

Finally, having studied the special cases of one and two 4-simplices, we now turn to the general case of a triangulation 
containing arbitrary many simplices. The proof of the existence of a solution to the simplicity constraint ([5]) is simple, 
and is performed by induction over the number of 4-simplices. As a base case, note that a single simplex can be 
labeled with 10 m’s and 10 L’s such that the simplicity constraint is satisfied, as already discussed above. To prove 
the inductive step, first assume that we have found a solution (m, L) that labels N 4-simplices. Then, assuming that 
one of those 4-simplices features a free boundary tetrahedron, we proceed to attach the {N + l)-th 4-simplex to that 
boundary tetrahedron. We then specify the labels (m, 1) of that additional 4-simplex to be equal to the corresponding 
labels of its neighbor, by the recipe discussed in the case of two 4-simplices above (namely (HU and (1181) ). The 
analysis of two 4-simplices above then guarantees that those labels also satisfy the additional six simplicity constraint 
equations, since the equations are the same for the two 4-simplices, and they are satisfied for one of them by induction 
hypothesis. Therefore, we have constructed an extended set of labels (m, L) which satisfies the simplicity constraint 
for -I- 1 4-simplices. This completes the proof of the induction step. 

At this point a few comments are in order. First, numerical search for solutions distinct from dnl), (HI has failed 
to find any, suggesting that this solution of the simplicity constraint equations is the only possible one, at least in 
the generic case. Second, as we can see from the construction, once we choose some particular labeling for the first 
4-simplex, all other 4-simplices in the triangulation are labeled in the same way, i.e. the triangulation consists of 
4-simplices which are all mutually identical. This has rather nontrivial physical consequences, because the state sum 
m reduces to the form 


zr™= i: i: i: 


JSR{Lc{m)) 


(19) 


where a counts the number of possible solutions for L’s, given a specific set of m’s, in the system m- From the 
state sum one can see that there are only ten degrees of freedom in the whole spacetime, specified by the identical 
labeling of all 4-simplices in T. This is unsatisfactory from the physical point of view, because we expect that the 
theory gives general relativity in the classical limit, and general relativity has two physical degrees of freedom (i.e. two 
graviton polarizations) per every point in space, which is certainly more than ten in total. This problem is commonly 


9 


resolved by passing to the formalism of “second quantization”, i.e. by taking an additional sum over various possible 
triangulations: 


= E E-- 

TeT TeT miez mioSZ a 

The set T is some nonempty set of inequivalent triangulations keeping the topology of the manifold the 

same. One can in principle attempt to discuss “all” possible triangulations, but it is often more useful to restrict 
to a certain class, according to various physical and mathematical criteria. Obviously, one of these criteria is the 
convergence of the state sum ([20|) , assuming of course that ([T^ is finite to begin with. 

We will not discuss the convergence issues or the choice of T in this paper. We just note that since T can be 
chosen to be suitably large, it can provide an arbitrarily large number of degrees of freedom in the theory, through 
different configurations of the triangulation. It is however important to emphasize that this is a highly nontrivial step, 
and there are a priori no guarantees that summing over some set of triangulations will indeed introduce the needed 
degrees of freedom into the theory. It may still happen that the model contains anomalies and the space of solutions 
is overconstrained. The presence of anomalies is an important open problem that should be studied in the future. 

As a final comment, note that the sums over m’s and a are labeling all 4-simplices at the same time, independently 
of the choice of T. Under the implicit assumption of uniform convergence of the sum over triangulations, one is 
therefore allowed to switch the order of summations in to obtain 

E ■■■ E EfE , (21) 

TTIiqGZ Ct \T’G7” / 

where the term in parentheses becomes very similar to the state sum discussed in the models of CDT approach to 
quantum gravity. In the next section we turn to the analysis of this relationship. 


E E 


JSR{La(m)) 


( 20 ) 


= E Zrp ^ 


V. RELATION TO CDT 


In the Causal Dynamical Triangulations approach (see 0 for a review), one constructs the state sum as 

^CDT = ^ , (22) 

TeT 

where the sum goes over a class of triangulations T specified by causality requirements, while all 4-simplices in T are 
isosceles, i.e. labeled by edge lengths a and 5, such that one can distinguish foliations of T into spacelike hypersurfaces 
labeled exclusively by edge lengths a, while each two hypersurfaces are connected by edge lengths b, which is usually 
chosen to be timelike rather than spacelike (although this is not mandatory). This means that all tetrahedra within 
a given hypersurface are equilateral, while all 4-simplices filling up a slice of spacetime between two hypersurfaces are 
isosceles. In particular, of all possible isosceles 4-simplices (for given a, b there are in total 40 inequivalent ones up to 
reflections and rotations, see Appendix [B]) , exactly two are used, depicted by the following two diagrams: 




In both diagrams, thick lines are of length a, while thin lines are of length b. In the notation of @ these two 4-simplices 
are denoted as (4,1) and (3, 2) respectively. 

The result of the previous section was the state sum (l?T|) . which is readily comparable to (l2^ . The main difference, 
however, lies in the fact that for the generic choice of m’s in (l?T]) all 4-simplices in the triangulation must be identical, 
as dictated by the simplicity constraint. This is in contrast to (l2^ . where two types of 4-simplices are being used. 
However, it turns out that this case is also covered by (EID, as we will now show. 

As was discussed in the previous section, one can label the whole triangulation by labeling one arbitrary initial 
4-simplex, and then employing the simplicity constraint to fix the labels for all other 4-simplices. In a generic case. 
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this procedure is unique. However, if the labeling of the initial 4-simplex is chosen in a special way, there can be more 
freedom for the choice of other 4-simplices. As an example, choose the labels of the initial 4-simplex such that it is 
isosceles, for example the left one of the two diagrams above. If one assigns length a to thick edges and length b to 
thin edges, the simplicity constraint for that 4-simplex reduces to the following two equations: 

7/2 |mi| = A^(a,a,a) = |m 2 | = A//(a, 6 ,&) = ^v^4&2 - . (23) 

This is because only triangles of type (a, a, a) and (a, 6 , b) appear in the 4-simplex. Moreover, it is obvious that the 
simplicity constraint can always be uniquely solved for a and &, given any choice of integers mi and m 2 : 


a = 2lp 




But now note that there is another 4-simplex that is also made only of those two triangles, namely the right 4- 
simplex in the diagram above. Its simplicity constraint is identical to (I23I1 . and thus already satisfied. Therefore, in 
this case the 4-simplices in the triangulation need not be all identical, which leaves more flexibility in the possible 
labeling of the triangulation. In this case, restricting m to isosceles 4-simplices of the above type, we obtain 

= E E (e 

mi£Zm2GZ VTgT / 


where the expression in the parentheses is now precisely the CDT state sum Restricting oneself to this choice 

of labeling, after “freezing out” the values for mi and m 2 , the spincube model with the strongly imposed simplicity 
constraint reduces to the CDT model of quantum gravity. 

It is worth noting that there are other isosceles 4-simplices for which the analogous construction may apply. In 
particular, there are two more pairs of 4-simplices which are labeled by only two edge lengths a, b and consist of 
only two types of triangles (see Appendix [B] for details). One can choose those alternative geometries to satisfy the 
simplicity constraint in a similar way, and obtain a state sum similar to CDT. This is however not equivalent to 
CDT, because such 4-simplices do not induce a natural foliation of spacetime into space and time, which is one of the 
important aspects of CDT. 

The presence of the CDT state sum (l22l) as one piece of the spincube state sum (l?T]) essentially means that the 
spincube model contains the CDT model as a special case. Specifically, if one studies only the isosceles configurations 
of 4-simplices within the spincube model, one will recover all the wealth of results that can be obtained within the 
CDT approach to quantum gravity. In addition the spincube model allows one to study non-CDT-like isosceles 
configurations mentioned above, as well as non-isosceles configurations, all of which can potentially give rise to novel 
effects, not present in the CDT approach. Extending the quantum gravity model from the CDT state sum to the 
spincube state sum is therefore useful, both as a way to study more general spacetime configurations, and as a way 
to study the relationship between the CDT formalism and the spinfoam formalism, which is a very interesting open 
problem. 


VI. TOPOLOGICAL RESTRICTIONS 

So far we have not discussed the topological and combinatorial properties of the manifold AI 4 and its triangulation 
There are two main questions to be addressed: 

(a) Does the form of the state sum (I^Tl) place any restrictions on the choice of the topology of the manifold A 44 ? 

(b) Does the choice of the topology of AI 4 place any restrictions on the form of (1^ ? 

Before we even begin discussing these questions, two remarks are necessary. First, it should be noted that we 
are interested in the realistic choice of four spacetime dimensions, which means that we are discussing 4-dimensional 
topological manifolds. The full classification of 4-dimensional manifolds is known to be an undecidable problem in the 
sense of Godel’s first incompleteness theorem SHI, and even among simply-connected 4-manifolds there are those 
which do not admit a triangulation to begin with, like the Eg manifold . It is thus clear that any conclusive 

analysis of the above two questions is hopeless. This section is therefore devoted to discussing some examples and 
reformulations of the questions, without giving any general answers. 
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Second, one should keep in mind that the topology of is kept fixed when summing over different triangulations 
in m- Given one initial triangulation T{M4), one can construct other, different triangulations of the manifold with 
the same topology using, for example, Pachner moves. Some suitable subset of these triangulations will define the 
domain T for the sum over triangulations. Additionally summing over different topologies of would represent the 
formalism of third quantization. 


topologies 


and is out of the scope of this paper. 

After these introductory remarks, we can concentrate on what can be said regarding the questions (a) and (b) 
above. To this end, it is instructive to rewrite the state sum m as 

= ZiQ + • • • + Zi, 

where 

a VtgT / 

Zg = ^ E f E , 

rhi,...,fh3 a \TeT / 



Here we have rearranged the summation over the ten ffi integers such that we can explicitly distinguish the pieces of 
the state sum where two, three, etc. of them are equal. The Ziq term corresponds to the generic situation when all 
m’s are mutually different. The Zi term corresponds to the case where all of m’s are equal. The Z 2 term contains 
the generalized CDT state sum (1241) . In relation to questions (a) and (b) above, it will turn out that the Ziq term 
will be most useful to discuss (b), while the Zi term is the most useful to discuss (a). 

The second ingredient we need is the 1-complex dual to the triangulation. It is constructed as follows: to each 
4-simplex in the triangulation one assigns a vertex, and to each tetrahedron common to two 4-simplices one assigns 
a link connecting two vertices. If the tetrahedron is on the boundary, then it belongs to only one 4-simplex, and 
the corresponding link is open-ended. Thus the 1-complex dual to the triangulation is some 5-valent graph with 
links encoding the adjacency relation among 4-simplices, while open-ended links encode boundary tetrahedra of the 
triangulation. 

Looking first at the generic case of Ziq, it was emphasized in section HVl that all 4-simplices must be identical, while 
having 5 different irregular tetrahedra on their boundaries. The 4-simplices are then glued to each other along the 
identical tetrahedra, keeping the simplicity constraint satisfied. However, there exist a possibility that, after a certain 
number of gluings, two distant parts of the triangulation are supposed to “meet” due to topological structure, and 
at the meeting point they may fail to have compatible tetrahedra for consistent gluing. In this way topology may 
invalidate the generic choice of m’s in Zio, and thus place a restriction the form of the state sum Z^'^ (see question 

(b)). ... 

This problem is actually equivalent to the problem of consistently coloring a 1-complex dual to the triangulation 
such that each link is colored by a number from the set {1,..., 5}, and such that at every vertex no two colors are 
repeated. We do not have a general answer to this question. Indeed, there may be 1-complex graphs for which such 
a coloring is possible, and graphs for which it is impossible. One obvious example of the former is the construction 
where we start from one vertex, connect five new vertices to its five links, and then keep connecting new vertices to 
new free links. But we do this in such a way that we never mutually connect two already existing free links, i.e. we 
never make a loop. The 1-complex graph constructed in such a way can always be colored in the required way — one 
starts by arbitrarily coloring the links of the initial vertex, and expands the coloring from there, keeping the coloring 
of links entering into each new vertex consistent with the previous choices. For the opposite example, we do not have 
any obvious examples, but it is certainly possible that for some graphs consistent coloring may be impossible. 
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The above analysis suggests that the topology of the manifold might exclude the generic piece Zio from the total 
state sum, but it might also accommodate it without problems. Regarding the question (b) above, this is the best 
possible answer one can give at this stage — depending on the actual choice of the topology of A^ 4 , there might or 
might not be some restrictions on the state sum. 

On the other extreme, let us look at Z\. This piece of the state sum contains equilateral 4-simplices. In an 
equilateral 4-simplex, there is only one type of tetrahedron available, and any two 4-simplices can be glued along any 
of their tetrahedra. In the language of the dual 1-complex, this means that we are supposed to color the links of 
the graph with a single color (vertices are now such that all “five” colors must be mutually identical). It is pretty 
obvious that every graph can be colored with a single color, without restrictions. This means that a triangulation 
corresponding to arbitrary topology can always be labeled with equilateral 4-simplices. 

The problem with this, however, is that equilateral 4-simplices can only be embedded into a spacetime of Euclidean 
signature, while we are interested in the realistic Lorentzian case. So the equilateral 4-simplices actually have to be 
excluded from the state sum on the grounds of desired spacetime signature. So the only available 4-simplices in Z\ 
are isosceles, with at least three different types of distinct tetrahedra. In the language of the dual 1-complex, this 
means that we are supposed to label the links of the graph with three distinct colors {1,2,3}, along with certain 
combinatorial rules for the vertices. Whether or not this is possible for a graph corresponding to the triangulation of 
an arbitrary topology, is again a problem for which we have no answer. Like in the generic case, one can certainly 
construct graphs such that this coloring is possible. In particular, every graph that can be consistently colored with 
five different colors (corresponding to the Z\q piece of the state sum), can also be colored using only three distinct 
colors (simply by identifying the fourth and fifth color with one of the previous three). For the opposite case we again 
do not have an explicit example, but it might exist. 

The analysis of the Zi case then suggests an answer to the question (a), albeit an inconclusive one — the state sum 
may place restrictions on the possible topologies of M 4 , but also maybe does not, depending on the existence of the 
dual 1 -complex which cannot be colored with three distinct colors. 

Finally, all remaining pieces of the state sum, Zg to Zg are “in between” the two extreme cases Zig and Zi, and 
add no new insight. As one goes from Zi to Zig, each piece Zk for higher k places more stringent rules for coloring 
the graph, ranging from three colors with repetitions at vertices to five colors without repetitions. In this sense their 
analysis does not add any novel information about possible topologies that is not already present in cases Zi and Zig. 

The inconclusiveness of the answers to questions (a) and (b) given by the qualitative analysis of this section 
might seem underwhelming and disappointing. One should therefore keep in mind that the undecidability of the 
classification of topological 4-manifolds severely limits what one can say about any general question involving the 
interaction between the structure of the state sum and topology, in the following way. Regarding the question (b), 
topology of AI 4 places restrictions on the state sum if no triangulation of M ,4 can accommodate the Zxg piece of the 
state sum. If we denote as S the set of all triangulations T that have some chosen fixed topology of AI 4 , 

S' = {T I T is homeomorphic to AI 4 }, 


then one can rewrite the statement that topology places restrictions on the state sum as follows: 


VT T G S ^ 


dual graph of T cannot be colored 
in a way compatible with Zig 


But this statement is impossible to prove, since it is undecidable whether any given triangulation has the topology of 
A 44 , so the set of triangulations with appropriate topology, S, cannot be effectively quantified over. In other words, 
the antecedent in the implication is undecidable, which makes the whole sentence undecidable. 

Regarding the question (a), the state sum places restrictions on the topology of AI 4 if there is no triangulation of 
AI 4 that can accommodate the Z\ piece of the state sum: 


VT T as ^ 


dual graph of T cannot be colored 
in a way compatible with Z\ 


Again, this is impossible to prove, for the same reason as before — it is undecidable whether or not any given 
triangulation has the topology of M. 4 ^ so the set of candidate triangulations cannot be quantified over. 

Nevertheless, the analysis presented above gives us some practical tools and techniques to study particular cases. 
For example, given some particular triangulation, one can always construct its dual I-complex, and determine whether 
or not it can be colored according to the rules given above. This question is decidable for any finite graph, since there 
are finitely many possible choices for the color of each link. One can simply try out all possible combinations of colors 
for all available links in the graph, and the algorithm is guaranteed to halt after a finite number of steps, since there 
are finite number of possibilities to test. 
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VII. CONCLUSIONS 

In this paper we have studied the possibility of strongly imposing the simplicity constraint in the spincube model of 
quantum gravity. After a short overview of the model, we have introduced two ways to impose the simplicity constraint, 
weakly and strongly. Weak constraint holds only in the classical limit of the theory, and naturally leads to the Regge 
model of quantum gravity. However, imposing the constraint strongly, such that it holds at the quantum level, leads 
to a nontrivial issue of consistency which depends on solving a complicated system of Diophantine equations. As the 
main result of the paper, it was shown that this Diophantine equations indeed do have solutions, and two classes of 
solutions were constructed. The first class represents the case when all 4-simplices in the triangulation are identical 
(up to reflection and rotation symmetries). The second class represents the case of isosceles 4-simplices, whose edge 
lengths take only two different values. Moreover, one can choose the labeling of those isosceles 4-simplices to be 
such that the theory is equivalent to the state sum of the CDT approach to quantum gravity. Finally, we have 
discussed the topological properties of the state sum, where the conclusions are very limited by the undecidability of 
the classification of topological 4-dimensional manifolds. 

We close with a few remarks. Note that the spincube model state sums (I^Tl) and (IMl) are more general than the 
CDT state sum ()22l) . in the sense that one sums over all possible choices for the triangle labels mi and m 2 , which are 
kept constant in the CDT model. The reason for this is that we have extended the spincube model to the formalism of 
second quantization — summing both over the field variables and the triangulations, rather than summing over only 
one of those. As was discussed, the second quantization formalism is necessary if one strongly imposes the simplicity 
constraint, since otherwise the model has too few degrees of freedom. Of course, passing on to the formalism of second 
quantization does not automatically guarantee the introduction of the correct number of degrees of freedom, since 
the theory may contain anomalies and the space of solutions may be overconstrained. This is common for both the 
spincube model and the CDT approach, and is an open problem that should be studied in the future. 

On the other hand, summing over both the field variables and the triangulations provides us with a better un¬ 
derstanding of the theory. Namely, one can resolve the long-standing problem of the relationship between the CDT 
approach and various covariant LQG approaches (spinfoam and spincube models), within the same unifying frame¬ 
work. Also, the formalism includes other types of solutions to the simplicity constraint and therefore other specific 
models inequivalent to CDT. In this sense, one can study those alternative models and compare them to the CDT 
approach. This opens many novel avenues for further research. 

Finally, one can study some nontrivial issues related to the classical limit. Namely, in CDT the classical limit 
is defined statistically, by averaging over many possible triangulations with a given (large) number of 4-simplices. 
Depending on the parameters in the action, one discovers three different phases in the theory, one of which can be 
related to the classical general relativity. In this setup the edge lengths of the 4-simplices are mutually equal and are 
kept fixed. On the other hand, in the spinfoam and spincube models, the classical limit is obtained via the asymptotic 
limit in which edge lengths, areas and volumes become very large, while the triangulation itself is kept fixed. It is an 
open question whether these two limiting procedures are equivalent or not, what is the meaning of each, and which 
procedure is more appropriate for the study of the classical theory emerging from the model. The second quantization 
formalism of the spincube model discussed in this paper provides one with the tools to study both limits within the 
same model. In particular, it would be interesting to reproduce the phase diagram of CDT geometries for a model 
consisting of irregular 4-simplices, where more than three edges (up to all ten) have different values. One could study 
the properties of resulting phases, and how they relate to the CDT phases as one varies the values of the various edge 
length parameters. Also, one could discuss whether in the limit of large edge lengths one always ends up in the same 
CDT phase or not, and why. All these questions are interesting to study, and will be addressed in future work. 
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Appendix A: State Sum Models 

A general state sum model in D dimensions is constructed as follows. Let T{A4d) be a triangulation of the manifold 
Md- Triangulation T contains vertices v, edges e, triangles A, tetrahedra r, and so on up to D-simplices a. Each 
of these objects is labeled with “colors” cj), a set of quantities that describe the fundamental variables of the model. 
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In addition, each type of objects in T is assigned an “amplitude” A, which enters the description of the dynamics of 
variables (j). One then writes the state sum Z as 

^=E n n • n 

{4>} v&T eeT creT 

The sum goes over the full domain of each color (j) living on each vertex, edge, triangle, tetrahedron or a higher order 
simplex. 

The above expression should be understood as a rigorous definition of a path integral 

Z = y Dcj) exp , (A2) 

where the D + 1 amplitude functions Av, Ae, .4a, • ■ •, Aa- partly enter into the definition of the path integral measure 
V(/) and partly into the definition of the action Roughly speaking, specifying the set of variables (j) and their 

domain is equivalent to the specification of the fields present in the model, while specifying the amplitudes A is 
equivalent to choosing their dynamics. 

The equation represents the discretization of (IA 2 I) in the sense of Feynman. In particular, if one constructs 
a state sum model over a one-dimensional manifold = K. (called the time axis) and chooses amplitudes Av and 
Ae conveniently, equation (lAll) reduces precisely to the Feynman’s textbook “definition by discretization” of the 
configuration-space path integral in quantum mechanics. Generalization of the same idea to multiple-dimensional 
manifold Ai d leads to (EU in the general case. 


Appendix B: Isosceles 4-simplices 

We define a given 4 -simplex to be called “isosceles” iff all its edges have length a or b, where by convention 
a,b & R'*’- All triangles in such a simplex are either isosceles or equilateral. Also by convention, we exclude fully 
equilateral 4 -simplices, by requiring that edges of both lengths a and b must be present in the 4 -simplex. 

It is an interesting question to determine all possible inequivalent isosceles 4 -simplices. There are in total 20 different 
ones, up to the a O & exchange symmetry and permutations of the vertices. If we label the vertices of a 4 -simplex as 
per the diagram below. 



1 2 


and if we order the ten edge lengths into a 10-tuple as {I12, hs, lu, ^15, hs, hi, hb, hi, hb, hb), the distinct labelings for 
the 20 isosceles 4 -simplices are given as: 


1 

(a. 

a. 

a, 

a. 

a. 

a, 

a, a, 

a,b ), 

11 

(0, 

0, 

a. 

b,a, 

o, 6,6,0,6), 

2 

(a. 

a. 

a, 

a. 

a 

a 

a, a 

b,b), 

12 

(0, 

a. 

a. 

6,6, 

6,0,6,0, a ), 

3 

(a. 

a, 

a, 

a. 

a 

a 

b, b, 

CL, u) , 

13 

(o. 

a. 

a. 

b,a, 

6,0,6,6, o), 

4 

(a. 

a 

a 

a. 

a 

a 

6, a 

b, b) , 

14 

(o. 

a. 

a. 

a, a 

a, 6,6,6,6), 

5 

(a. 

a 

a 

b. 

a. 

a 

b,b, 

a, a ), 

15 

(0 

a 

b. 

6,6, 

o, 6,6,0, a) 

6 

(a. 

a 

a 

a. 

a 

a 

b, b, 

0,6), 

16 

(0 

a 

a 

6,0 

b, 6,6,6, a) 

7 

(a 

a 

a 

a. 

a 

a 

a, b 

6,6), 

17 

(0 

a 

a 

6,0 

a, 6,6,6,6), 

8 

(6, 

b. 

a, 

a. 

a. 

K 

a, a, 

a, a ), 

18 

(0 

a 

a 

b,a 

6,0,6,6,6), 

9 

(0 

a 

a 

b. 

a 

a 

b, a, 

6,6), 

19 

(0 

a 

a 

a, a 

6,6,6,6,6), 

10 

(a 

a 

a 

b. 

a 

a 

b, b, 

0,6), 

20 

(0 

a 

a 

6,6, 

6,0,6,0,6). 


The above table has been computer-generated by brute force counting of all possibilities. If we mark all a-edges with 
thin lines and all 6-edges with thick lines, these 4 -simplices can be drawn as follows, respectively: 
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All the above 4-siinplices have different 4-volumes, tetrahedra volumes and triangle areas. To these one should also 
add 20 more 4-simplices, which are obtained by exchanging 006 , i.e. by switching thick and thin lines. This gives 
a total of 40 inequivalent isosceles 4-simplices. 

There are four different types of triangles that make up isosceles 4-simplices. These are (a, a, a), (a,a,b), (a,b,b) 
and (b, b, b) triangles. Most of the 4-simplices contain either three or sometimes all four types of triangles, making 
them unsuitable for satisfying the simplicity constraint, as explained in the main text. Nevertheless, there are five 
pairs of 4-simplices which are made up only of two types of triangles. These are: 



containing only (a, a, b) and (b, b, b), 



containing only (a, a, a) and (a, a, b), 



, containing only (a, a, b) and (a, b, b), 


and their a b duals (the third pair is self-dual). Of these pairs, one can use the first pair (or their dual) to construct a 
triangulation which features 3-dimensional hypersurfaces made of equilateral tetrahedra, such that these hypersurfaces 
are separated by “equal distances” everywhere. Thus one can establish the triangulation-induced foliation of spacetime 
into space and time. This property is exploited by the CDT approach to quantum gravity. The remaining pairs above 
lack this property, because they mix edges of size a and 6 in a way that does not admit a nice 3-1-1 foliation. 


Appendix C: Asymptotics of the weak-constrained Regge action 


Here we give the proof of the statement from the main text that the action (HU becomes equal to the proper Regge 
action in the asymptotic limit of large edge lengths. As noted below equation (HU, the limit is defined by rescaling 
all edge lengths If, by a common factor k, and then expanding the resulting action into the asymptotic series when 
k ^ 00 . Therefore, we start from HU evaluated at klf, 


SR{kl) = 


E 

AeT 


—I^Auiklf^A) 


S^ikl), 


and expand it into the asymptotic series as follows. First we use the fact that for any x £ R, the floor function can 
be written in the form 


[xj = X — i?(x), 

where the remainder i?(x) is always such that R{x) £ [0,1). Taking x to be AH{kl)/jlp and applying the above 
formula to our action, the first term can be recognized as the proper Regge action ©, while the second term contains 
the reminder. We thus obtain 


Snikl) = 5Regge(fc0 - ^ R{kl)6A{kl ), 
AgT 


(Cl) 


where 

R{kl) = R (^^Anikl)^ , 

and therefore 


0 < R{kl) < 1 


(C2) 
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for every € R. 

Next we need to recall from geometry that the deficit angles are homogeneous functions of order zero in the edge 
lengths, while Heron formula for triangle area is a homogeneous function of order two in the edge lengths, 

dA{kl) = Sa^I) , AH{kl) = eAH{l). 

Therefore, the Regge action ([7]) scales as SB.egge{kl) = k'^S^eggeil), so we can rewrite (|Cip as 


— 'S^Regge(^0 


1 - 




Due to (IC2I1 . the second term on the right-hand side scales as 0{l/k^) in the limit k —)> oo. Therefore, in this limit 
the action (ED does indeed become asymptotically equal to the Regge action ©, 


SRikl) = Sneggeikl) 




{k —>■ oo), 


as stated in the main text. 


Appendix D: Numerical analysis of simplicity constraint solutions 

In section IIVI we discussed the system of simplicity constraint equations (fT51) for a single 4-simplex, and claimed 
that the set of solutions is quite rich, based on the numerical Monte-Carlo analysis. Here we provide some details of 
that analysis. 

We begin by noting that, for every triangle A, its simplicity constraint equation ([SD can be rewritten in the form 

|”^a| = —i^Ah{Ii,12 , h) = Ah (Li, L2, L 3 ) , 

where Lk = h/^lp G R are the dimensionless edge lengths of the triangle A. This is due to the fact that the Heron 
formula for the triangle area is a homogeneous function of order two in the edge lengths. Therefore, if we label the 
vertices of a 4-simplex as in the diagram 



we can write the simplicity constraint system of equations explicitly as follows: 


|to123| 
|to124| 
|to125| 
|to134| 
1^135 I 
|toi45| 
|to234| 

1 ^ 235 ! 

1^2451 

1 ^ 345 ! 


Ah {Li 2 , L 13) L 23 ) 
Ah {Li 2 , Lii, L 24 ) 
Ah {Li 2 , Li5 , L 25 ) 
Ah {Li3, -bi4, 1 ^ 34 ) 
Ah (Lis, Lib, L 35 ) 

Ah {Li4, Li 5 , L 45 ) 

Ah {L 23 , L 24 , L 34 ) 

Ah {L 23 , L 25 , L 35 ) 

Ah {L 24 , L 25 , L 45 ) 

Ah {L 34 , L 35 , L 45 ) 


(Dl) 


Here, each triangle and each edge is labeled with the vertices that belong to it. 

The algorithm we used for the analysis of the solutions of the system dniD goes as follows. First, we randomly choose 
the ten integers ttia S {!,..., 50}, where restricting to positive integers does not lead to loss of generality. Given 
this data, we proceed to look for all numerical solutions of m for the ten edge lengths Lg G [0,100]. By randomly 
seeding 1000 initial guesses for ten edges Lg from the given domain, we employ the implementation of the Newton’s 
method to find solutions of dnn, up to machine precision. All found solutions are memorized and counted. For each 
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solution, we calculate the Euclidean volume of the corresponding 4-simplex, using the Cayley-Menger determinant 
formula 



0 

1 

1 

1 

1 

1 


1 

0 

^12 

/■2 

-^13 

/■2 

-^14 


1 

1 

/■2 

^12 

0 

^23 

^24 

^23 

9216 

1 

^13 

/■2 

-^23 

0 

^34 

^33 


1 


-^24 

/■2 

-^34 

0 

-^45 


1 

■'^15 

-^25 

^35 

-^45 

0 


(D2) 


Only solutions having different 4-volumes are retained. After the seed of 1000 initial edge lengths has been exhausted, 
the algorithm returns to the beginning, randomly choosing another set of integers ttia , and repeats the calculation. 

One should note that the Newton’s method in general does not guarantee to hnd all possible solutions in a given 
range, since some of the solutions may only be reached by choosing a seed outside the given domain. Nevertheless, 
for our purposes it is not necessary to find all solutions, but merely some of them, since our goal is to demonstrate 
that the system m generically has many solutions. 

The results of the analysis are as follows. We have counted the number M of randomly chosen integers ttia that have 
n different solutions for the ten edge lengths L^, such that the 4-volume of each solution is different from the others, 
and such that the 4-simplex can be embedded into a 4-dimensional space with Lorentz signature. The latter condition 
is achieved if the Cayley-Menger determinant (imi) is negative. The function M (n) thus gives us the distribution 
of choices of tua over the number of distinct Lorentzian solutions for Lg of the simplicity constraint system dnu. 
Normalized over 1000 choices of toa, the obtained distribution is: 


n 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

M{n) 

0 

5 

15 

37 

51 

87 

106 

137 

124 

112 

107 

89 

64 

28 

16 

8 

7 

4 

2 

1 

0 


This can be represented as a histogram, in the following way: 



As we can see, for a random choice of ten ttia, the simplicity constraint system jnu is most likely to have somewhere 
between 6 and 12 solutions with Lorentzian signature. No choices of have been found to have only one, or more 
than 20 solutions. The record number of 20 solutions have been found for the following choice: 


TO123 = 28, 
TO124 = 42 , 
TO125 = 43 , 
TOi34 = 32 , 
TO135 = 26, 


mi45 = 27, 
W234 = 42 , 
W235 = 38 , 

^245 = 28 , 
771345 = 24 . 


(D3) 


The corresponding solutions for as well as the values of the Ej, are as follows (note that these results have been 
obtained with 64-bit machine precision, but are quoted here only up to three decimal places, for brevity): 
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a 

Li2 

.^13 

Li4 


L23 

to 

L25 

.^^34 

L35 

L45 


1 

9.188 

6.866 

9.647 

15.263 

8.572 

15.298 

9.790 

10.301 

9.907 

7.216 

-96358.047 

2 

9.751 

6.180 

15.423 

9.195 

9.420 

9.135 

15.179 

11.610 

8.769 

7.755 

-83095.409 

3 

9.841 

6.254 

15.930 

10.998 

9.213 

9.267 

9.294 

11.839 

15.816 

6.433 

-138374.305 

4 

9.936 

6.799 

9.492 

10.013 

8.329 

10.153 

9.948 

12.365 

15.309 

5.825 

-14626.743 

5 

9.995 

6.191 

10.364 

9.425 

9.251 

9.278 

10.576 

11.682 

8.615 

18.949 

-29435.832 

6 

10.010 

6.334 

10.261 

15.153 

8.992 

9.341 

8.940 

12.964 

10.411 

6.579 

-65535.329 

7 

10.012 

6.292 

10.173 

15.024 

9.063 

9.409 

8.896 

11.887 

10.365 

6.576 

-43342.444 

8 

10.152 

7.352 

8.705 

9.266 

7.644 

15.291 

10.615 

11.388 

15.097 

6.465 

-70393.956 

9 

10.362 

18.623 

9.328 

10.066 

9.212 

9.937 

9.522 

10.528 

9.378 

18.492 

-428991.593 

10 

10.669 

18.040 

8.539 

8.688 

8.425 

10.780 

16.077 

10.868 

10.257 

6.831 

-184052.701 

11 

10.809 

6.742 

9.996 

15.572 

15.987 

8.991 

8.363 

10.156 

10.273 

7.101 

-121442.024 

12 

10.888 

17.907 

8.255 

8.851 

8.106 

15.831 

10.485 

11.059 

9.978 

6.932 

-139620.742 

13 

11.257 

5.519 

16.580 

10.493 

10.168 

8.263 

8.653 

13.213 

9.466 

7.365 

-7145.581 

14 

11.451 

6.478 

10.104 

8.715 

16.443 

8.608 

10.288 

10.800 

8.600 

17.808 

-162923.501 

15 

11.518 

6.032 

15.321 

9.311 

15.853 

7.551 

9.551 

11.737 

8.979 

7.573 

-67437.306 

16 

11.710 

7.059 

9.385 

7.464 

8.084 

19.156 

12.714 

13.118 

9.413 

7.386 

-2106.96 

17 

11.875 

19.723 

8.718 

9.650 

8.671 

9.798 

9.092 

12.076 

10.768 

6.333 

-26923.715 

18 

13.967 

6.434 

10.134 

8.083 

19.418 

8.362 

10.683 

13.002 

10.254 

6.697 

-4065.232 

19 

14.069 

7.238 

8.847 

9.425 

8.957 

9.618 

9.214 

11.595 

15.150 

6.327 

-68517.683 

20 

14.739 

6.577 

9.842 

8.854 

10.111 

8.807 

9.966 

10.986 

8.320 

17.665 

-119422.745 


One can see that all solutions describe 4-simplices with different 4-volumes. Moreover, it is straightforward to verify, 
using Heron’s formula, that for each choice of edge lengths, all ten triangle areas are integers (lD3l) . 

Finally, we should note that the only purpose of the results obtained in this analysis is to illustrate that the set of 
solutions of the system m is indeed quite rich, as claimed in the main text. 
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